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Abstract 

The evolution of a hadronic system after its chemical decomposition is described 
through a model that conserves the hadronic multiplicities to their values at chemical 
freeze-out. In the partition function describing the model all known hadronic reso- 
nances with masses up to 2400 MeV have been included. The state of the system is 
found as function of temperature and the corresponding baryon density is evaluated. 
The baryon density at thermal decoupling is also computed. 
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1. Introduction 

Thermal approaches have extensively been used to describe the particle multiplicities 
which emerge from high energy collisions [1-16]. The results of such approaches are satisfac- 
tory since they are able to predict quite accurately a large number of different experimentally 
measured hadronic multiplicities as function of a few thermodynamic variables, such as tem- 
perature, volume and chemical potentials. 

The extracted parameters from such approaches define the "chemical freeze-out" , i.e. the 
point where the chemical composition of the system that produces the particle multiplicities 
is fixed. Another point called "thermal freeze-out" can also be defined. This second point 
is associated with the particle momentum distribution which is measured experimentally. 
After this point this distribution remains fixed and the particles no longer interact among 
themselves. 

There is evidence that in a lot of circumstances the two points corresponding to the same 
system do not occur at the same temperature. Generally the freeze-out temperature is lower 
than the chemical freeze-out one. Since the particles are measured once they have reached 
the experimental apparatus, after any kind of interaction among themselves has ceased and 
since the thermodynamic parameters of the chemical freeze-out point predict quite well all 
those particles, one has to infer that all these abundances will have to remain fixed through 
the whole process between chemical and thermal freeze-out. 

In this paper the main focus will be to construct a model of relativistic hadronic parti- 
cles formulated in the grand canonical ensemble that will be able to conserve the particle 
multiplicities between points which correspond to different temperatures. Of course the two 
points with greatest interest are the chemical and the freeze-out point. The question that 
arises then is why there is need for another model and why not use one of the existing ther- 
mal models. In these models the particles are described as thermally equilibrated entities. 
In some of these models the hadrons are non-interacting particles [2-7] and in others a kind 
of interaction among them has been included [1,8-16]. But usually the free thermodynamic 
parameters are the volume V, the temperature T and a set of a few chemical potentials 
each of which is associated with the conservation of quantum numbers like baryon number 
B, charge Q, strangeness S, etc. During the evolution of the system among states with 
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different temperatures its content can alter. This is done by adjusting the chemical poten- 
tials appropriately so as to keep the relevant quantum numbers fixed. The only parameter 
that remains unfixed after the conservation of quantum numbers is the volume. If there is 
need to fix the particle numbers as well it is clear that this cannot be accomplished with 
the existing parameters. Even if it is assumed that conservation of all the particle numbers 
automatically conserves the quantum numbers, meaning that only the conservation of the 
particles is enough, this limits the application of the above models to only situations where 
the number of the particle entities is equal to the number of the chemical potentials plus 
oneS. But generally the hadrons which have to be considered are more numerous than the 
quantum numbers potentials. 

As it is evident these thermal models cannot evolve the system from chemical to thermal 
freeze-out. In this work the necessity to have fixed particle numbers will be used to construct 
a model which will determine the evolution of the hadronic system after its chemical freeze- 
out. This newly constructed model will coincide at chemical freeze-out with "Ideal Hadron 
Gas" model (IHG) [2-5], one of the aforementioned thermal models, which is formulated in 
the grand canonical ensemble and describes hadrons as relativistic non-interacting particles. 

2. Formulation of the model 

Before discussing the new model, the IHG model will be presented. In the context of 
IHG the grand canonical partition function, formulated in the Boltzmann approximation, 
has the form 

In Z (V, T, {\})ihg =VJ2 X QNi Z H tj ( T ) = V^QN^^SiMM^-) > « 

i j i j Z7r 1 

where i runs over all hadronic families such as mesons, N Baryons, A Baryons, etc. and 
j represents the specific member of the family with degeneracy factor g>jj and mass m%j. 
\QNi stands for the product of all the fugacities associated with the particular family. These 
fugacities can either be quantum numbers fugacities related to Baryon number, Strangeness, 
etc. or to quark flavour!!. 

1 In such a case the quantum number potentials and the volume would have to be adjusted to keep the 

particles fixed between points with different temperatures. 

2 For example, for S~ Baryons, Xqn would read ^b^q 1 ^ 2 ! 2 or ^d^ll 2 - One can look in [16], eq. 14, to 



2 



One can evaluate particle abundances if the above partition function is extended by the 
introduction of a fugacity Ajj for every particle. After calculating the particle number one 
has to set in IHG Ay = 1 [17], so again the particle number is only expressed as function of 
the quantum numbers fugacities. 

Now, if someone wishes to keep the particle numbers fixed it is only natural to extend 
(1) by the use of fugacities Ay (corresponding to particle numbers), but with the difference 
that the constraint Ay = 1 will not be imposed. This model will be called Fixed Particle 
Numbers (FPN) model and accordingly its partition function will depend on Ay's 

T m ■■ 

\nZ(V,T,{X}) FPN = Vj2\ Hij Z Hij (T) ee VE^, 3 ^<^(f) , (2) 

ij ij 

where A# 4j is product of quantum numbers as well as particle number fugacities!. Since the 
experimentally measured multiplicities usually contain feeding from the decay of resonances, 
all known hadrons with masses up to 2400 MeV have been included in the FPN partition 
function. The same hadrons have also been included to IHG partition function (1). The 
mean particle number can be evaluated through the relation 

d\nZ(V,T,{X}) FPN 



< Nij >FPN— Ay- 



dXij 



(3) 



where {A} ^ Ay means that for the evaluation of the partial derivative all fugacities except 
Ay are considered as constants. 

The next point that has to be elucidated is at what values the particle numbers will stay 
fixed. Since after chemical freeze-out these values do not alter it is useful to formulate FPN 
so as to keep the particle numbers fixed at their chemical freeze-out values. One has to 
remember that the chemical freeze-out values are extracted from a thermal model, like IHG, 
by a successful fit to the experimentally measured values. The particle numbers can then 
be fixed at the values calculated through IHG for the chemical freeze-out thermodynamic 
variables. Thus the results of IHG and FPN should coincide at chemical freeze-out point. 



find out how the two sets of fugacities are related. 

3 For example, for S(1530)~ Baryon, Xh would read AsAg 1 Ag 2 7g A S ( 1530 )- 



3 



So it has to be required! 

< N {j >ihg=< >fpn^ VX QNi Z Hij {T) = V'\ Hij Z H ..{T') & 

Hij ~ V'Z Hij (T>) ' 1 } 

The above equation can be used to calculate the total product of fugacities consisting 
of quantum numbers fugacities and of hadron fugacities at temperature T'. It has to be 
pointed out that it is not possible to evaluate each quantum number fugacity separately, but 
this is irrelevant since the full product of fugacities can be calculated. Another focal point is 
that all quantum numbers are automatically conserved as linear combination of the particle 
numbers. 

In the right hand side of eq. (4) the only term which is left undetermined after the impo- 
sition of the conservation of particle numbers is the multiplicand factor yr. So an additional 
constraint has to be applied. For example conservation of entropy can be assumed!. 

The entropy of the system can in general be calculated from! 

s = - ( a ±IM^IAm) , (5) 



dT J 



where \i represents the chemical potential associated with fugacity A = exp(yu/T). Applying 
(5) to the partition function (2)! the constraint of fixed entropy will read 

S = S' & 

# \nZ IHG {V,T,{n}) + VT^QN^P- - VT^X QNi ^Z Hij (T) = 

ij ij 

= In Z FPN (V\ T', {//}) + V'T' £ X HiJ 0Zh ^ T) - V'T' E A^. ^Z Hij (T') . 



The primed variables in this paper will generally be related to subsequent points of the chemical freeze- 
out point. 

5 A lot of authors assume isentropic evolution of the system, e.g. see [18]. 

6 The symbol of entropy is tilded in order not to be confused with the symbol of Strangeness. K can be 
set equal to one. 

7 Since the IHG and FPN partition functions coincide at chemical freeze-out point, the IHG partition 
function (1), where A^- = 1, can be used for the evaluation of the entropy at this point. 
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With the use of (4) the last equation becomes 
VT E \ QNi dZH ^ T) - V E X QNi H*QN>)Z Hij (T) = 



dT 

f V'Z Hij (T>) " dT' ' ' ^ V'Z Hij (T>) { V'Z Hij (T>) ] 

ij ij 

-1 2^ A Q^z 7^77 gj7, 2^ ^QNi^Hij [J- ) m( — — — ) . 

(6) 



Setting x = (6) can be solved for a; to give 



exp 



l^ij ^QNi^H i:j [J- ) ln[ <Z H (T') / + J ^ij df 1 ^ij A Q^i z H (T 7 ) 9T> 



(7) 



Equation (7) can be used to evaluate the volume expansion ratio as the system has cooled 
to a temperature T' less than the chemical freeze-out temperature T. With the use of the 
same equation, quantities like the baryon density of the system can be calculated at T' . 
The baryon chemical potential at which the system is found at temperature T' cannot be 
calculated separately from the rest of chemical potentials in the context of FPN. But the 
baryon density has no problem to be evaluated. One has to remember that baryon number 
is also fixed with the imposition of the constraints (4). So 

< B' > _ < B > ch _ V < B> ch n c g /oA 

" li FPN ~ 



V V V V x FPN ' 

where rf^ is the baryon density calculated at chemical freeze-out. 

3. Application 

The newly constructed model, FPN, can describe thermally equilibrated hadronic systems 
with fixed particle numbers when their temperature is known. As an example the systems 
formed at different interactions at SPS will be considered. For this reason the chemical freeze- 
out parameters obtained for these systems through fits to their experimentally measured 



values will be used. These parameters are listed in Table 1 along with the references where 
they can be found. From a variety of thermal analyses performed by different authors the 
particular ones have been chosen because they allow for partial strangeness equilibrium 
(7 a ^ 1) and they use most recent available values for the experimentally measured hadronic 
multiplicities. The values of Table 1 are then taken, for each interaction separately, as input 
to the equations < S >= and 2 <q> = to determine the rest of the fugacities. Thus the 
whole set of chemical freeze-out parameters (T, fj, B , /j,q, fis,7s) are calculated and also the 
products of fugacities Aqjv, in (1) are also set. 

Giving different values to temperature T, equation (8) can be used to calculate the 
corresponding baryon density. The resulting paths for FPN for S + S, S + Ag and Pb + Pb 
interactions are shown in Figure 1 with solid curves. For the Pb + Pb interaction the thermal 
freeze-out temperature is calculated in Refs. [20] and [21]. For these values baryon density at 
thermal freeze-out n l g er can be evaluated. The results are listed in the last column of Table 
2. The path for Pb + Pb is followed until the lower temperature (of the two given in Refs. 
[20,21]) is reached. The points which correspond to the thermal freeze-out temperatures of 
these references are depicted with squares on the FPN curve. 

FPN has the unique attribute to conserve each particle species separately. So it is not 
possible to compare its results directly with another thermal model. In order to have a 
general view we shall depict on the graphs with the FPN results IHG states for different 
temperatures. The IHG model to be used only conserves the quantum numbers < B >, 
< Q > and < S > and the entropy < S >. The IHG states are represented by a dotted 
curve in Figure 1. It has to be pointed out that IHG does not take the system from chemical 
to thermal freeze-out. The system remains in the context of IHG all the time at chemical 
equilibrium where the particles can transform into one another. 

The usefulness of Figure 1 is that, as can be observed, a "soup" of fixed particles cools 
more efficiently than a chemically equilibrated IHG state. This means that at the same 
baryon density (which through the conservation of the baryon number is equivalent to equal 
volumes) smaller temperature corresponds to the FPN state than to the IHG state. The 
reason is that IHG only conserves a few quantum numbers. As temperature drops the 

8 /3 is fixed from the baryon number and charge of the participant nucleons, e.g. see [16]. 
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number of particles diminishes but without affecting the preservation of quantum numbers. 
For example an equal reduction to the number of protons and antiprotons will not affect the 
conservation of baryon number. In the FPN case, on the contrary, no particle number is 
allowed to diminish. So in the IHG state the available energy has to be distributed among 
less particles than FPN, and their mean kinetic energy has to be greater, leading to greater 
temperature. 

In Figure 2 the FPN ratios x = V'/V, where V is the chemical freeze-out volume, 
are plotted as function of temperature for the three SPS interactions. For comparison the 
ratios x = VI hg /Vjhg f° r the particular IHG model discussed above are plotted with dotted 
curves. The two volumes in the IHG ratios correspond to chemically equilibrated states at 
different temperatures. V' lHG and Vjhg are calculated for the same temperatures as V and 
V respectively. 

Also, in order to show the effect of hadrons with large masses, FPN calculations with 
hadrons only up to the Delta mass have been included in figures 1 and 2. The chemical 
freeze-out parameters for Pb + Pb of Table 1 have been used for these calculations. From 
figure 1 it is evident that the baryon density calculated through the truncated version of 
FPN is considerably less than the baryon density calculated at the same temperature with 
the FPN model using all the hadrons. The expansion ratio x for the truncated FPN, though, 
is close to the calculated ratio with FPN using all the hadrons, as it can be seen from figure 
2. 

4. Conclusion 

After chemical freeze-out the collisions among hadrons that compose the hadronic gas 
can no longer change its chemical composition. Following this requirement a non-interacting 
hadron gas model (FPN) has been presented that keeps the multiplicity of every particle 
fixed to the value dictated by the chemical freeze-out conditions. In the context of FPN 
the constraints of conservation of quantum numbers are broken up to a larger number of 
constraints, these of conservation of particle numbers. The chemical potentials of quantum 
numbers are no longer "good" variables to describe the evolution of the system. Of course 
the fugacities of particle numbers used as variables in FPN are not "free" parameters. Their 
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values are fixed from the given set of the quantum numbers fugacities at chemical freeze-out. 
So the evolution of a hadronic system is described as function of temperature (after imposing 
conservation of entropy). This is done for three SPS interactions. 

Following this evolution and using values of thermal freeze-out temperature extracted for 
the Pb + Pb interaction the baryon density at freeze-out is evaluated. As the temperature at 
thermal decoupling for various interactions can be calculated using transverse mass spectra 
or HBT analysis [22] the same procedure can be applied to evaluate the corresponding baryon 
density before free streaming for these interactions. 

In this paper the IHG has been used as the thermal model which would coincide with FPN 
at chemical freeze-out point. Any other thermal model, interacting or non-interacting, can 
also be used in place of IHG and, with the use of particle fugacities, a model that conserves 
particle multiplicities can also be formed. Recently in model II of [23]! the particle fugacities 
have been used to conserve particle multiplicities. This a three dimensional hydrodynamic 
model describing radial and elliptic flow but it includes hadrons with masses only up to Delta 
mass and it is restricted to zero baryon chemical potential. 

Acknowledgement I would like to thank Professor N. G. Antoniou and Professor 
C. N. Ktorides for reviewing the manuscript and for useful remarks. 

References 

[1] J. Cleymans, H. Satz, Z. Phys. C 57 (1993) 135. 

[2] J. Cleymans, K. Redlich, H. Satz, E. Suhonen, Z. Phys. C 58 (1993) 347. 

[3] J. Sollfrank, M. Gazdzicki, U. Heinz, J. Rafelski, Z. Phys. C 61 (1994) 659. 

[4] J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank, J. Rafelski, Phys. Rev. D 51 (1995) 3408. 

[5] A. D. Panagiotou, G. Mavromanolakis, J. Tzoulis, Phys. Rev. C 53 (1996) 1353. 

[6] F. Becattini, U. Heinz, Z. Phys. C 76 (1997) 269. 
9 [23] appeared during the refereeing period of the present paper. 



8 



[7] F. Becattini, M. Gazdzicki, J. Sollfrank, Eur. Phys. J. C 5 (1998) 143. 

[8] D. H. Rischke, M. I. Gorenstein, H. Stocker, W. Greiner, Z. Phys. C 51 (1991) 485. 

[9] R. Hagedorn, Nuovo Cimento Suppl. Ill (1965) 147. 

[10] R. Hagedorn, J. Ranft, Nuovo Cimento Suppl. VI (1968) 169; R. Hagedorn, Nuovo 
Cimento Suppl. VI (1968) 311. 

[11] R. Hagedorn, Nuovo Cimento LVI A (1968) 1027. 

[12] R. Hagedorn, J. Rafelski, Phys. Lett. B 97 (1980) 136. 

[13] A. S. Kapoyannis, C. N. Ktorides, A. D. Panagiotou, J. Phys. G 23 (1997) 1921. 
[14] A. S. Kapoyannis, C. N. Ktorides, A. D. Panagiotou, Phys. Rev. D 58 (1998) 034009. 
[15] A. S. Kapoyannis, C. N. Ktorides, A. D. Panagiotou, Phys. Rev. C 58 (1998) 2879. 
[16] A. S. Kapoyannis, C. N. Ktorides, A. D. Panagiotou, Eur. Phys. J. C 14 (2000) 299. 
[17] R. Hagedorn, K. Redlich, Z. Phys. C 27 (1985) 541. 
[18] J. Cleymans, K. Redlich, Phys. Rev. C 60 (1999) 054908. 

[19] F. Becattini, J. Cleymans, A. Keranen, E.Suhonen, K. Redlich, Phys. Rev. C 64 (2001) 
024901. 

[20] H. Appelshauser et al, NA49 Collaboration, Nucl. Phys. A 638 (1998) 91c; Eur. Phys. 
J. C 2 (1998) 661. 



[21] J. R. Nix et al., |nucl-th/980l045 



[22] U. Heinz, J. Phys. G 25 (1999) 263; R. Stock, Nucl. Phys. A 661 (1999) 282c. 

[23] T. Hirano, K. Tsuda, talk given at 30th International Workshop on Gross Properties 
of Nuclei and Nuclear Excitation: Hirschegg 2002, in Hirschegg 2002, Ultrarelativistic 
heavy-ion collisions 152-157, e-Print Archive: |nucl-th /0202033 . 



9 



Figure Captions 

Figure 1 Contours (solid lines) that follow hadronic systems after chemical freeze-out on 
(T,n B ) plane for 3 interactions at SPS, calculated through FPN (model of Fixed Particle 
Numbers). On the same graph points at chemical equilibrium (dotted lines) calculated 
through an IHG model that conserves baryon number, charge, strangeness and entropy are 
depicted. The slashed curve represents calculations with FPN including hadrons with masses 
only up to the Delta mass (FPN (a)) and using as chemical freeze-out parameters the ones 
that correspond to Pb + Pb interaction of Table 1. 

Figure 2 The ratios of the volumes V of the SPS hadronic systems at a certain temperature 
to their volumes V ch at chemical freeze-out calculated for FPN (solid lines). The dotted lines 
correspond to the ratios of volumes of chemically equilibrated states of IHG model of Figure 
1 calculated at the same temperatures as the volumes of the FPN ratios. The slashed curve 
represents calculations with model FPN (a) of Figure 1. 
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Table Captions 

Table 1 Chemical freeze-out parameters calculated for different interactions at SPS and the 
corresponding references. 

Table 2 Thermal freeze-out temperature calculated in two different references for the Pb + 
Pb interaction and the corresponding computation of baryon density through FPN. The 
upper errors of baryon density correspond to the upper errors of temperature. The same is 
true for the lower errors. 



Experiment 


T ch (MeV) 


H c £(Mt 


iV) 


if 


Reference 


S+S 200 A ■ GeV 


180.5 ± 10.9 


220.2 ± 


18.0 


0.747 ±0.048 


[7,18] 


S+Ag 200 A ■ GeV 


178.9 ±8.1 


241.5 ± 


14.5 


0.711 ±0.063 


[7,18] 


Pb+Pb 158 A ■ GeV 


174.7 ±6.7 


240 ± 


14 


0.900 ± 0.049 


[19] 



Table 1. 



Experiment 


T ther (MeV) Reference 


riB er (M~ 3 ) 


Pb+Pb 158 A ■ GeV 
Pb+Pb 158 A ■ GeV 


120 ± 12 [20] 
95.8 ±3.5 [21] 


o.o99iffi 
o.o627i8:88I 



Table 2. 
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